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This  combined  Quarterly  Progress  Report  is  abbreviated.  The  research  done  in  this  period  which 
has  appeared  as  ERL  reports  have  been  circulated.  We  apologize  for  mailing  this  report  written  well 
after  the  work.  The  time  and  support  needed  to  produce  QPR’s  as  in  the  past  were  not  available.  Also, 
Prof.  Birdsall  was  Joint  Institute  of  Fusion  Theory  Visiting  Professor  at  Nagoya  University  during  most 
of  this  period. 

Our  QPR  history  dates  from  the  late  1960’s,  and  has  served  us  well  as  a  vehicle  to  communicate 
our  research  results  and  as  practice  in  writing  for  our  research  students,  prior  to  journal  article  and 
thesis  writing.  The  latter  has  helped  produce  an  average  of  3-5  publications  per  Ph  ■  D.  -.uuwi,  Ul 
•ime  of  graduation.  However,  times  change,  support  changes,  and  emphasis  changes. 

We  list  our  publications  for  all  of  1988,  1  Journal  article,  4  ERL  reports,  9  talks  and  5  invited 
talks.  Abstracts  of  the  IEEE  talks  and  APS/DPP  talks  are  provided. 

Reprints  of  the  Journal  article  and  the  ERL  reports  have  already  been  sent. 

C.K.  Birdsall 
Principal  Investigator 
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Publications  for  1988 


Journal  Articles 


Niels  F.  Otani  and  Bruce  I.  Cohen,  "Effect  of  Large-Amplitude  Perpendicularly  Propagating  Radio  Fre¬ 
quency  Waves  on  the  Interchange  Instability,"  Phys.  Fluids,  31,  pp.  158-176,  January  1988. 


ERL  Reports 


K.  Theilhaber,  "Vorte:  Formation  and  Particle  Transport  in  a  Cross-Field  Plasma  Sheath,"  University  of 
California,  Berkeley,  Memorandum  No.  UCB/ERL  M88/21,  March  20,  1988. 

Lou  Ann  Schwager  and  C.  K.  Birdsall,  "Collector  and  Source  Sheaths  of  a  Finite  Ion  Temperature  Plas¬ 
ma,"  University  of  California,  Berkeley,  Memorandum  No.  UCR/ERL  M8/23,  April  13,  1988. 

Lou  Ann  Schwager,  "Effects  of  Secondary  Electron  Emission  on  the  Collector  and  Source  Sheaths  of  a 
Finite  Ion  Temperature  Plasma,"  University  of  California,  Berkeley,  Memorandum  No.  UCB/ERL 
M88/24,  April  13,  1988. 

Lou  Ann  Schwager,  "Effects  of  Ion  Reflection  on  the  Collector  and  Source  Sheaths  of  a  Finite  Ion 
Temperature  Plasma,"  University  of  California,  Berkeley,  Memorandum  No.  UCB/ERL  M88/25,  April 
13,  1988. 


Talks,  Conference  Proceedings 

JCRF/Edge  Physics  Workshop.  March  30-April  1,  1988,  Boulder,  Colorado. 

W.  S.  Lawson  and  C.  K.  Birdsall,  "Undriven  Plasma  Wall  Interaction  Simulations,  Showing  Tur¬ 
bulence  with  and  without  an  Initial  Vacuum  Gap." 

W.  S.  Lawson  and  C.  K  Birdsall,  "Antenna  Driven  Plasma  Wall  Interaction  Simulation,  Showing 
Local  Turbulence  and  About  3  Times  Larger  Flux  to  the  Wall." 

IEEE  Conference  on  Plasma  Sciences.  June  6-8,  1988,  Seattle,  Washington. 

C.K.  Birdsall,  "Serendipity  is  no  accident,  even  in  plasma  research",  Plasma  Science  and  Applica¬ 
tions  Award  acceptance  address. 

W.  S.  Lawson,  M.  A.  Lieberman,  and  C.  K.  Birdsall,  "Electron  Dynamics  of  RF  Driven  Parallel 
Plane  Reactor." 

A.  Friedman,  S.  L.  Ray,  C.  K.  Birdsall,  and  S.  E.  Parker,  "Multi-Scale  Particle-in-Cell  Plasma 
Simulation:  Timestep  Control  Criteria  and  some  Tests." 

APS  Division  of  Plasma  Physics  annual  meeting,  Hollywood,  Florida,  October  31-Novembcr  4,  1988. 

S.E.  Parker,  "A  Proposed  Particle-In-Cell  Method  for  Modeling  Small  Angle  Coulomb  Collisions 
in  Plasmas." 

C.K.  Birdsall,  and  W.  S.  Lawson,  "Computer  simulation  of  Plasma  Response  Near  an  RF 
Antenna." 

IJ.  Morey,  R.W.  Boswell,  and  C.K.  Birdsall,  "Particle  Simulation  of  a  Low  Pressure  RF 
Discharge." 

P.G.  Gray,  W.S.  Lawson,  and  C.K.  Birdsall,  "Heat  Flow  Between  Species  in  a  One-Dimensional 
Simulation  Plasma." 
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Talk 

C.K.  Birdsall,  "Plasma-Sheath-Surface  Dynamics  via  Particle  Simulations,"  June  9-10,  1988,  at  NSF 
Workshop  on  New  Directions  in  Plasma  Engineering,  UC  Berkeley. 

Invited  Talks  (in  Japan) 

C.K.  Birdsall,  "Vortex  formation  and  particle  transport  in  a  cross-field  plasma  sheath,"  (Co-author  K. 
Thcilhaber),  September  26-27,  1988,  at  Plasma  Sheath  and  Potential  Formation  Seminar,  Sendai  Japan. 

C.K.  Birdsall,  "Plasma  sheath  small-signal  RF  impedances  as  obtained  from  simulations  of  a  planar 
device,”  October  28,  1988,  Seminar  at  Inst.  Plasma  Physics,  Nagoya  University,  Japan. 

C.K.  Birdsall,  "Source  and  collector  sheaths  in  a  bounded  plasma  device."  (Co-author  L.A.  Schwagei , 
and  "The  magnetized  plasma  sheath  interacting  self-consistently  with  an  absorbing  wall;  Kelvin- 
Helmholtz  instability  growth  with  saturation  as  a  dynamic  steady  state,  producing  Bohm  diffusion," 
(Co-author  K.  Theilhaber),  December  8,  1988,  Inst,  of  Plasma  Physics,  Nagoya  University,  Japan. 
(The  latter  also  presented  at  Hiroshima  University  a  week  later.) 

C.K.  Birdsall,  "Hands-on  demonstration  on  a  personal  computer  of  our  periodic  code  ESI  and  of  our 
bounded  code  PDW1,"  (Co-author  John  Verboncoeur)  December  9,  1988,  Inst,  of  Plasma  Physics, 
Nagoya  University,  Japan. 


APS  Division  of  Plasma  Physics  Annual  Meeting 
November,  1988 


1R  19  Computer  Simulation  of  Plasma  Response  Near  an  RF 
Antenna.  C.  K.  BERDSAXL  Electronics  Research  Lab,  U.  C.  Berkeley 
W.  S.  LAWSON  Courant  Institute,  New  York  University  —  A  two- 
dimensional  electrostatic  code  is  used  to  investigate  the  response  of 
a  magnetized  plasma  to  RF  drive  at  ion  cyclotron  frequencies  as  a 
function  of  drive  frequency  and  amplitude.  The  plasma  parameters 
are  chosen  to  model  a  Tokamak  edge  plasma  as  closely  as  possible. 
The  electrostatic  approximation  is  of  marginal  validity,  but  the  results 
should  be  indicative  of  the  important  effects  in  the  electromagnetic  case. 
The  plasma  characteristics  of  interest  are  the  ion  flux  to  the  wall,  and 
the  energy  of  the  ions  that  strike  the  wall.  Our  initial  results  show 
that  the  flux  to  the  wall  increases  with  the  drive  amplitude,  and  that 
the  energy  of  ions  striking  the  wall  is  a  strong  function  of  frequency. 
Simulations  at  a  drive  frequency  of  10/3  of  the  ion  cyclotron  frequency 
showed  no  increase  in  temperature,  but  those  at  5/3  showed  an  order 
of  magnitude  increase  or  more.  Further  results  will  be  presented. 

This  work  was  supported  by  TJSDoE  contract  DE-AC03-87ER53254 


2P  15  Heat  Flow  Between  Species  in  a  One-Dimensional  Sim¬ 
ulation  Plasma-  P.  G-  GRAY  Dartmouth  College,  W.  S.  LAWSON 
Courant  Institute,  NYU,  and  C.  K.  BIRDSALL  UC  Berkeley  —  We 
apply  the  kinetic  theory  of  a  one-dimensional  plasma  as  worked  out  by 
Eldridge  and  Feix  [Ref.  l]  to  the  problem  of  heat  flow  between  species  in 
a  simulation  plasma.  It  is  necessary  to  define  a  simulation  temperature 
which  may  be  very  different  from  the  intended  physical  temperature. 
Heat  flows  in  accordance  with  this  simulation  temperature,  allowing 
beat  to  flow  from  one  species  to  another  which  was  intended  to  be 
hotter.  This  effect  can  be  a  serious  pitfall  to  simulators  who  model  a 
minority  species  as  a  large  number  of  sub-particles  each  of  which  has  a 
fraction  of  the  charge  of  the  particles  of  the  other  species. 

1  O.  C.  Eldridge  and  M.  Feix,  Phys.  Fluids  8,  398  (1963). 


2Q  8  A  Proposed  Particle-In- Cell  Method  for  Modeling 
Small  Angle  Coulomb  Collisions  in  Plasmas.*  S.  E.  Parker,  Elec¬ 
tronics  Research  Laboratory,  U.  C.  Berkeley  —  We  are  developing  a 

"•'•’’•'••i  •luxlel  small  angle  Coulomb  collisions  seif-ronsistently  This 
mciitou  may  be  used  in  standard  PIC  plasma  simulations  to  include  col¬ 
lision-  .  ->r  as  .in  alternative  to  solving  the  Fokker-Planck  equation  using 
finite  ilifference  methods.  The  distribution  functions  for  both  ions  and 
electrons  are  represented  by  a  large  number  of  particles.  The  particle 
velocities  change  as  a  fnnctio.i  of  the  drag  force  and  the  diffusion  in  ve¬ 
locity,  which  is  represented  by  a  random  walk  process.  This  is  analogous 
to  previous  Monte-Carlo  methods1,  except  we  calculate  the  drag  force 
and  the  diffusion  tensor  self-consistently.  The  particles  are  weighted  to 
a  grid  in  velocity  space.  The  associated  “Poisson’s  equations"  are  soived 
for  the  Rosenbluth  potentials,  from  which  the  drag  force  and  diffusion 
tensor  .are  obtained.  This  poster  will  outline  the  proposed  method. 

Work  performed  for  USDoE  under  contract  FG03-86ER220. 

1  T.  D.  Ronglien  end  T.  A  Cutler,  .Vuclear  Fusion  20,  1003  (1080' 


2Q  18  Particle  Simulation  of  a  Low  Pressure  RF  Discharge 
LJ.  Morey,  B_W.  Boswell*  ,  C.K.  Birds aE,  Electronics  Research  Lab¬ 
oratory,  U.  C.  Berkeley  —  RF  discharges  have  been  simulated  with 
a  hounded,  one  dimensional  particle- in-cell  model  which  includes  both 
ionization  and  secondary  electron  emission.  One  electrode  is  grounded 
and  the  other  Is  made  to  oscillate  between  ±200V  with  a  frequency  of 
20MHz.  The  simulations  are  started  with  a  small  number  of  cold  ion 
and  electron  pain  and  run  until  equilibrium  is  reached.  At  equilibrium 
the  ions  and  electrons  react  to  the  average  and  instantaneous  poten¬ 
tials  respectively1  since  u,,  <  uiT;  <  uip,.  A  greater  proportion  of  the 
ionization  occurs  within  the  sheaths  at  higher  pressure,  but  there  is  an 
optimnmpressure  at  which  the  highest  plasma  densities  are  achieved. 
It  was  sis"  found  that  the  discharge  could  be  sustained  without  any 
secondary  electron  emission,  and  that  the  effects  of  secondary  electron 
emission  due  to  the  ions  decreases  with  pressure. 

*  st  the  Australian  National  University. 

1  R.W.  BoeweU  end  IJ.Morey,  Appl.  Pbrs-  Let t.  52,  21  (1988). 
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